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Abstract 

We define and evaluate on the lattice the amputated 3-gluon vertex function 
in momentum space. We give numerical results for 16 3 x 40 and 24 3 x 40 
quenched lattices at f3 = 6.0. A good numerical signal is obtained, at the 
price of enforcing the gauge-fixing condition with high accuracy. By comparing 
results from two different lattice volumes, we try to investigate the crucial issue 
of finite volume effects. We also outline a method for the lattice evaluation of 
the QCD running coupling constant as defined from the 3-gluon vertex, while 
being aware that a realistic calculation will require larger f3 values and very 
high statistics. 
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I. INTRODUCTION 



The task of measuring the QCD coupling constant at energies of the order of the b— quark 
mass is a major challenge for the lattice community. Many different methods have been 
proposed in this respect, based on the study of the interquark static potential, charmonium 
spectra and other quantites [1-3]. In spite of the success of some of these methods (at 
least when applied in the quenched approximation), it would be desirable to investigate a 
more fundamental, field-theoretic definition of the coupling constant, arising from the lattice 
study of some fundamental QCD interaction vertex, and compare the results with what is 
obtained from different methods. We discuss here some studies in this direction [4]. 

On general grounds, the lattice study of gluon (and quark) correlation functions is very 
interesting because it provides a valuable tool to gain insight into the non-perturbative QCD 
dynamics at the fundamental level. In fact, in this approach one deals with the fundamental 
degrees of freedom of the theory, assuming we are sufficiently close to the continuum limit, 
and the results can be compared in principle to continuum non-perturbative predictions 
from Schwinger-Dyson equations. 

Using a lattice definition of the gluon, one can evaluate any unrenormalized, complete 
ra-point Green's function (we omit Lorentz indices): 

fti,.--,^) = <Ai at (x 1 ),...,Ai at (x n )>, (1) 
or, in momentum space: 

G\: ) t (p 1 ,..., Pn )S\p 1 + ...+p n ). (2) 

Once G\ll(p 2 ) and G\al{pi,P2,P3) are determined, many interesting issues can be investi- 
gated. The first one is the non-perturbative behaviour of the gluon propagator 
which has been analyzed by many authors on the lattice [5-9], yet it is still not completely 

(2) (3) 

understood. Secondly, from G la j. and G la j. one may define non-perturbatively the lattice 
version of the amputated, 1PI 3-gluon vertex function: 
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rS2(pi,P2,P3) = g^p^ps) x n [gSCp?)]" 1 • (3) 



1=1 



Evaluating the above quantity is the crucial step in order to define a QCD coupling constant 
from the lattice 3-gluon vertex, as we explain in the following [4]. 

In order to determine a convenient kinematical setup, we first consider the general form 
of the continuum, off-shell 3-gluon vertex [10], which can be written as 

r (3) = -p(3) X r (3) L /.x 

cont ct{3^ cont ct{3^ ' cont aP^f \ J 

where 

r< fl^c t p 1 (Pi,P2,P3) = F(pl,pl;pl)(g al3 p x -p 2 -pipP2cc) B * 

+H [-g a p + ^ (pi 7 p 2a P3/3 ~ Pip P2 7 P3a)] + cyclic permutations, (5) 



wi 



ith 



#y = (Pi-y P2-P3- P2 1 Pi ■ Ps), (6) 



anc 



T< ?lt L c t p-y(Pi,P2,P3) = A(pI,p 2 2 ;p 2 3 ) g a p (pi ~P2) 1 + B(pl,p 2 2 ;pl) g aP ( Pl + p 2 ) 1 

+ C(p1,pI;pI) (Pip P2a - gap Pl ■ Pi) (Pl - P2) 1 + ~ (Pl-y P2a) P3P + Pip P2 1 P3a) + C.p. (7) 

In the above expressions the scalar functions F, A and C are symmetric in their first two 
arguments, B is antisymmetric in the first two arguments, S is antisymmetric under exchange 
of any pair of arguments and H is totally symmetric in p\-,p\-,p\- 

(3) 

Summarizing, T cont a p J contains 6 independent scalar functions, but for the purpose of 
computing the coupling constant renormalization one only needs to determine the function 
which multiplies the tree-level vertex, namely A(pl : pl;pl). Of course this is the only one 
which is divergent when the UV cutoff is removed. 

If one evaluates T^J nt a ^ at the asymmetric point defined by 

" = 7^, Pi = -Ps, P2 = 0, p 3 = pp, (8) 



then it can be written as 

r ( fJ ntaPa (- Pp ,0, Pp ) = 2A( P 2 ) Pp . (9) 

The above expression is proportional to the continuum tree-level vertex evaluated at (8). 

One can show that the leading term of the 1-loop lattice calculation of the 3-gluon vertex 
for the same kinematics is indeed consistent with (9) [11], thus when calculating the lattice 
vertex function T la j. at the point (8) one can set (neglecting terms of higher order in the 
external momentum): 

r£ aPa {-Pp^Pp) = 2 A lat (p 2 ,a) pp, (10) 

where a is the lattice spacing. Then, when a — > 0, one has 

Aiat{p 2 ,a)\ p 2 =q 2 = Z~ 1 (a 2 q 2 ) g (a), (11) 

for a generic momentum q 2 (g (a) is the bare lattice coupling constant). Then we obtain 
the relation 

z;\a 2 q 2 ) 9o (a) = -±- rg a p a (-Pe,0,Pe)\v>=*- (12) 

Analogously, one defines the gluon wavefunction renormalization constant Za from the re- 
lation 

G\H ^(P 2 )\ p i= q i = TM Z A (a 2 q 2 ) 1 (13) 

with T^y^q) being the projector on transverse fields (we will work in the Landau gauge). 
Finally, following Ref. [11], we define the renormalized, "running" coupling as 

g R (q 2 ) = Z 3 A /2 (a 2 q 2 ) Z~\a 2 q 2 ) g (a). (14) 

At this point, a standard 1-loop calculation may be performed to relate the coupling constant 
as defined above to more popular continuum schemes, like for example MS, etc. 

One drawback of the above procedure is that one needs an independent determination 
of the value of the inverse lattice step size ^ in units of energy. It has been recently pointed 
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out that a very precise determination of such a parameter can now be obtained from the 
analysis of the charmonium spectrum on the lattice [2]. Also, the determination of ^ from 
string tension measurements has now reached a high level of accuracy [12]. 

On the other hand, a very attractive feature of the above method is that in principle it 
can be applied to the unquenched case without any modification. 

It is important to emphasize that in the above procedure one is assuming that lattice 
correlation functions in momentum space can be effectively parametrized according to con- 
tinuum formulas. In this respect, provided that (3 is sufficiently large, it is crucial to study 
the role of finite lattice spacing and finite volume artifacts. The former ones should be 
under control as long as the momenta that we inject lie well below the value of the inverse 
lattice spacing. In other words, when injecting the momentum k x = 2 L V ™ along the lattice 
x direction, with n an integer, we require k x to satisfy the condition k x <C \- This gives 

(15) 

Finite volume effects are potentially more dangerous since the kinematics we have chosen 
for the definition of the vertex requires the evaluation of the lattice Fourier transform at 
zero momentum. In this exploratory work we will focus on the investigation of finite volume 
effects. 

As far as the actual calculation of the coupling constant is concerned, one is interested in 
evaluating it at momentum scales well above 1 GeV, where a comparison to the perturbative 
result may be attempted and confinement effects can be neglected. In particular, evaluating 
the coupling at q = mj ~ 5 GeV would allow a comparison to the experimental result [13]. 

On the lattices which are presently available to us, one cannot inject momenta above 
1.5 GeV and still assume that finite lattice spacing artifacts are negligible. For this rea- 
son, the results presented here are not suitable for the extraction of the coupling constant. 
Nonetheless, they provide evidence that the 3-gluon vertex function can be effectively mea- 
sured so that our method for the evaluation of the running coupling may be applicable at 
larger (3. 



II. THE GLUON PROPAGATOR 

The first step in our program is the evaluation of G lat (p 2 ). The lattice gluon field can 
be defined as [5]: 

2iag 3 \ 2iag J 

For the purpose of defining the amputated vertex function, we consider a set of 25 configu- 
rations on a 16 3 X 40 lattice at (3 = 6.0 (the same considered in [8]) and 33 configurations 
on a 24 3 X 40 at the same value of (3 (an enlargement of a set considered in Ref. [8], where 
a better gauge-fixing accuracy has been achieved). 

By comparing the results from such sets we hope to gain insight into the role of finite 
volume effects. 

All the calculations are performed after gauge-fixing to the so-called minimal Landau 
gauge [14,15] (see Ref. [8] for a short review), implemented through the iterative minimiza- 
tion of 

Hu[g] = "TtE R e Tr {^{n) + U?(n - A)) , (17) 

where V is the lattice volume and U 9 indicates the gauge-transformed link U 9 (n) = 
g{n)U fj,{n)g\n + fi). The typical precision that we reach for such a minimization is 10~ 5 . 
In Fig. 1 we show our results for the momentum space gluon propagator on the 16 3 X 40 
lattice. 



III. THE 3-GLUON VERTEX 

After calculating the propagator, we proceed to the evaluation of the complete 3-point 
function G la [ a f3 a {—pp,0,pp). As our main purpose is to test the quality of the signal and 
the occurrence of finite volume effects, we proceed as follows. First we set a = i, i = 
1, . . . , 3, (3 = 4, injecting momentum in the longer (time) lattice direction, and we call p t the 
momentum. In this case we can safely inject a few units of momentum while satisfying the 
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requirement given by Eq. (15), but finite volume effects are expected to be most severe as 
the zero-momentum Fourier transform is performed along all the (shorter) spatial axes. To 
improve statistics one observes that from Eq. (7) it follows that at the asymmetric point: 

r Slt 1 4 l (-Pt,0,p t ) = T^t 2 4 2 (-Pt,0,p t ) 

= T { l t 3 4 3 (-Pt,0,Pt). (18) 

Obviously an analogous relation holds for G^ 3 \ so that we can average G^ and I^ 3 ) over 
the spatial Lorentz components. 

After this calculation, we change the kinematics. We define an index which runs on 
all lattice directions but the i one. Then we set a = /3 = z, z = 1, . . . , 3 and we inject 
a momentum pi in the i direction. In this case not only we expect finite volume effects to 
be less severe, but we also have better statistics, as the symmetries of (7) allow to average 
over 9 distinct vertex functions. Finally, we compare the results from the 2 sets of lattices. 

Consider first the 16 3 X 40 lattice. The complete 3-point function G\ a l(—p,0,p) vs. p 
is shown in Fig. 2 for both the kinematical arrangements described above. The points 
corresponding to space and time momenta are plotted with different symbols. It appears 
that for large values of the momentum p, G^ 3 \— p, 0,p) gets damped, due to the effect of 
the propagators lying on the external legs. We observe that the data points are roughly 
consistent with an odd function of the injected momentum, as expected. In this respect, the 
points obtained from spatial momentum, besides having smaller error bars than the others, 
turn out to be antisymmetric to a better accuracy. 

In order to define the amputated vertex, we multiply the complete 3-point function shown 
above by the inverse propagators, according to Eq. (3). The resulting function T la j. (—p, 0,p) 
is shown in Fig. 3, where the error is a jackknife one. 

As the vertex must be an odd function of the momentum, we can further improve statis- 
tics by antisymmetrizing. The results are shown in Fig. 4. 

Turning now to the bigger lattice, i.e. the 24 3 X 40 one, we show in Fig. 5 the amputated 
vertex function after antisymmetrization. 
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By comparing the results in Figs. 4 and 5 it appears that, at the level of precision that we 
can presently reach, finite volume effects are not detected, as the data points corresponding 
to the injection of momentum in time are consistent between the two lattices. 

As expected, the data points corresponding to spatial momenta are the most stable. We 
collect such points from both lattices in Fig. 6. 

We observe that, in spite of the better statistics, the results from the bigger lattice are 
noisier. This may be related to an insufficient gauge-fixing accuracy. To investigate this issue 
we have performed some tests on a single configuration, studying the stability of our results 
vs. additional gauge-fixing. Such tests seem to indicate that the gauge-fixing accuracy was 
sufficient for our purposes, although this point deserves further investigation. 

In this connection, it is worth remarking that the accuracy of the numerical gauge-fixing 
is a very important parameter in our calculations. In fact, it turns out that for the purpose 
of evaluating gluon 3-point functions one needs a higher accuracy than when calculating 
the propagator. To illustrate this point, in Figs. 7 and 8 we plot the propagator and the 
complete 3-point function respectively, as evaluated on a single configuration of the 24 3 X 40 
lattice, for different levels of gauge-fixing precision. 

From Fig. 7 one can see that the propagator from 600 to 1400 gauge-fixing sweeps is 
virtually unchanged. On the other hand, Fig. 8 shows that for the 3-point function 600 
gauge-fixing sweeps are definitely not enough. 

IV. CONCLUSIONS 

We have discussed a method to measure on the lattice the 3-gluon vertex function in 
momentum space and to extract from it the running coupling constant. We have given 
numerical results at (3 = 6.0 for the vertex function. Such results are characterized by a 
good signal-to-noise ratio and show the symmetries expected. As a consequence, it appears 
that the operation of "amputating" at the non-perturbative level the lattice 3-point function 
in momentum space is feasible and leads to sensible numerical results, even with limited 
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statistics. 

By comparing results for two different lattice sizes it appears that finite volume effects 
do not play a significant role at our current level of precision. 

As we already remarked, our results are not yet suitable for a realistic determination 
of the QCD running coupling constant, both because of the level of statistical noise and 
because of the momentum range under consideration. Nonetheless, they suggest that such 
a task may be feasible in a high statistics calculation at weaker couplings. 
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FIGURES 

FIG. 1. G\ 2 J t (p 2 ) vs. p in GeV on the 16 3 x 40 lattice at /3 = 6.0. We assume a" 1 = 2.1 GeV. 

FIG. 2. (— p, 0,p) vs. p in GeV on the 16 3 x 40 lattice. Crosses and diamonds correspond 

to the injection of momentum in space and time directions, respectively. 

FIG. 3. Ij j (— p, 0,p) vs. p in Gev on the 16 3 x 40 lattice. Crosses and diamonds correspond 
to the injection of momentum in space and time directions, respectively. 

FIG. 4. Same as in Fig. 3, after imposing antisymmetry. 



FIG. 5. Same as in Fig. 4, for the 24 3 x 40 lattice. 

(3) 

FIG. 6. T y la j. (—p,0,p) vs. p in Gev for spatial momenta. Crosses and diamonds correspond to 
the 24 3 x 40 and the 16 3 x 40 lattices, respectively. 

FIG. 7. G^ t (p 2 ) vs. p in GeV calculated on a single configuration of the 24 3 x 40 lattice with 
different gauge-fixing accuracies. 

FIG. 8. G^ t (—p, 0,p) vs. p in GeV calculated on a single configuration of the 24 3 x 40 lattice 
with different gauge-fixing accuracies. 
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